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ABSTRACT 

Context. NGC 4441 is a candidate for a merger between a spiral and an elliptical galaxy (S+E merger), because it shows typical tidal 

structures such as an optical tail and two shells. With a far-infrared luminosity of ~ 5 • 10' Lq this galaxy belongs to the class of 

moderate luminosity mergers, in which the merging process induces (if at all) only a moderate starburst. 

Aims. The study of the atomic gas content allows us to investigate the merger history and the impact on the star formation. In particular, 

it is not clear from simulations whether an S+E merger leads to a gas concentration in the nucleus, resulting in a starburst, or whether 

the gas is spread out and therefore too diffuse for new star forming regions. 

Mefhods. We used the Westerbork Radio Synthesis Telescope to observe the properties of the Hi. By using this interferometer, we 

are able to study the large-scale H i distribution and kinematics with high spatial and velocity resolution. 

Results. We found two Hi tails extending out to more than 40kpc. In a central disc, the gas shows a fairly regular rotation pattern 

indicating that the gas started to settle after the merger. The total Hi mass adds up to 1.5 • 10' Mq. By comparing the high resolution 

Hi maps with deep optical images, differences between the stellar and gaseous tidal features are apparent, which indicates an S+E 

merger origin. 
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1. Introduction 

Extreme mergers between two large disc galaxies, leading to a 
super-starburst and to ULIRGs (Sanders & Mirabel 1996) have 
been at the focus of research, but those galaxies are rare. Small 
spirals and ellipticals dominate the galaxy population in many 
environments, but mergers between them, so-called S+E merg- 
ers, have been poorly studied so far Models of S+E merg- 
ers are in strong disagreement concerning the prediction of en- 
hanc ed star formation induced by the merger. IWeil & Hernquisj 
(1199 3) used a smoothed particle hydrodynamics approach (SPH) 
to model the gas. With this fluid description, the stellar and 
gaseous components are quickly segregated. While the stars 
form features like shells, the gas congregates in the centre of 
the remnant galaxy, leading to a strong gas concentration and 
thus resulting in a s t arburs t, similar to ULIRGs. Simulations by 
iKoiima & Noguchil (1 1997b however predict a dispersion of the 
gas clouds which might not lead to a starburst at all, because the 
density of the gas is too low to collapse and form new stars. They 
model the gas component as a system of inelastic cloud particles 
dissipating kinetic energy in mutual collisions. If the timescale 
for cloud-cloud collisions is relatively long, the cloud system be- 
haves like a collisionless system, thus like the stellar component. 
Therefore, the gas and stellar distribution are similar. While the 
contrast of the structures is somewhat higher for the stellar com- 
ponent, both show shell-like features, as were also found for the 
stars in the models of Weil & Hernquist ( 1993j)- 
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Simulations by'Naab & Burkert' ('2000^,'Naab & Burked (1200 ih . 
Khochfar & Burkert (2003) and Naab et al. (2006) revealed the 
importance of S+E mergers fo r the formation of bright ellip- 
ticals. As .Khochfar & Burkera (l2003h claim, these aspects are, 
however, poorly studied so far, and detailed observational inves- 
tigations are needed besides more numerical simulations to bet- 
ter understand the process of those mergers. 
One 'prototypical' S+E merger candidate is NGC 4194, the 
Medusa. In the optical we see a diffuse tail going to the north and 
on the opposite side 2 shells are visible, as predicted by models. 
Molecular gas is found in a continuous distribution out to 4.7 kpc 
from the cent re, i.e much more spread out than in the cas e 
of a ULIRG (lAalto & Htittemeisteii l2000t lAalto et al.1 l2001h . 
However, this galaxy is clearly undergoi ng an intense st arburst 
phase, albeit not as strong as in ULIRGs dWeistrop et al.l l2004). 
Here we present high resolution interferometric H i data of an- 
other S+E merger candidate, NGC 4441. This galaxy is mor- 
phologically very similar to the Medusa. It possesses one op- 
tical tidal tail and two bluish shells on the opposite side, which 
are, however, brighter, i.e. more evolved, than t he shells in the 
Medusa. The main body has an elliptical shape ('Bergvalll ll981h 
with a small dust layer throu gh the centre on the minor axis. 
Based on spectroscopic data iBergvalll (Il98lh claims that this 
galaxy is a merger where most of the gas has been already used 
for star formation, which means only little ongoing star forma- 
tion. However, there must have been a period of enhanced star 
formation, since the stell ar population i s younger than that of a 
normal elliptical galaxy (lBergvaIllll98lh . Our own optical spec- 
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Table 1. Basic properties of NGC 4441 . The distance is based on 
Ho = 75kms-'Mpc-'. 



property 


RA (2000) 


12:27:20.3 


DEC (2000) 


+64:48:05 


VopUei (kms"') 


2722 


D (Mpc) 


36 


type 


SAB0+ pec 


Lb(IO'Lo) 


10.1 


LpiRdO'Le) 


5.4 


1' 


10.5 kpc 



Table 2. Parameters for the WSRT observations. 



observing parameters 

date 

centre frequency 

total bandwidth 

number of channels 

velocity resolution (km s"') 

rms noise level (robust=0.5) (mjy beam"') 

beamwidth (") 

primary beam 

primary calibrator 



28.2.2003 

1408.9 

20 MHz 

1024 

4.12 

0.4 

20x19 

0.6° 

3C48 



tra confirm this an d we estimate that a moderate starburst oc- 
cured ~ 1 Gyr ago JManthev et al.ll2005l) . 

The observations presented here were obtained with the 
Westerbork Synthesis Radio Telescope. We refer to the analysis 
of optical long-slit spectra for metalicity estimations, but these 
data are presented in detail in Manthev et al. (2008c). 
In Table[T]general information on NGC 4441 is given. In Section 
2 we give basic information about the observations obtained, in 
Section 3 we describe the data reduction for the interferometric 
data. Results are presented in Section 4, including derived prop- 
erties like the H i gas mass and the star formation rate based on 
the 20 cm continuum flux. In Section 5 we discuss the obser- 
vations with particular emphasis on models of mergers between 
Ellipticals and Spirals and we compare NGC 4441 with the S+E 
merger prototype NGC 4194 (the Medusa). Finally, in Section 6 
we summarise the observations and results. 



2. Observations 

2.1. WSRT 

The Hi observations with the Westerbork Synthesis Radio 
Telescope (WSRT) were carried out in February 2003. Table|2] 
summarises the most important observing parameters. We ob- 
served NGC 4441 in one 12 h run in the maxi-short array con- 
figuration. This is best suitable for observations of an extended 
source within a single 12 hours track. The quality of the data 
is good, no interference problems occurred. For bandpass and 
flux calibration, the calibrator 3C48 was observed before the 
observations. Note that no secondary phase calibrator was ob- 
served. Instead, phase calibration at the WSRT is usually done 
using self-calibration. The chosen setup with a total bandwidth 
of 20 MHz (corresponding to ~4000km s"') and a velocity res- 
olution of 4km s"' channel"' was capable to detect both narrow 
features in velocity space and a broad total velocity range due to 
tidal disruption. With the chosen setup, two polarisations were 
observed, with 1024 channels available for each of them. 



2.2. Optical imaging 

We present here an optical R band image which was taken with 
CAFOS at the 2.2 m telescope on Calar Alto, Spain (CAHA) 
in April 2004. To study very faint morphological features, we 
obtained deep imaging with an integration time of 60 minutes. 
Standard reduction using IRAF was applied to the data. Here we 
use the optical image to compare the stellar and gaseous dis- 
tribution. For further analysis of the optical data, we refer to 
[Manthe y et al.. (2008c.) . 

3. Data Reduction 

The Hi data were reduced using MIRIAEQ] with some addi- 
tional tasks written by T. Oosterloo which are specific to the 
Westerbork interferometefl After inspecting and, where neces- 
sary, flagging the calibrators, the WSRT data have to be cali- 
brated first for changes in the system temperature. Since no sec- 
ondary calibrator was observed, after flux and bandpass calibra- 
tion a continuum image was created to perform self-calibration 
on continuum sources within the field. To achieve a proper cal- 
ibration, it is necessary to deeply CLEAN the continuum im- 
age. In an iterative process, first a mask was created which 
included only real emission to define the CLEANing compo- 
nents. Second, the image was CLEANed using the mask and 
subsequently a new, deeper mask could be generated for a bet- 
ter CLEANing. Once the image was finally CLEANed, self- 
calibration was performed, and the loop started again with cre- 
ating a mask, CLEANing and self-calibration. The data qual- 
ity was high enough that only 3 iterations were necessary for 
a proper calibration. The results were then applied to the tar- 
get galaxy. After continuum subtraction, a datacube was created 
to which we applied Hanning smoothing, leading to a veloc- 
ity resolution of 8 km s"'. The datacube was convolved with a 
Gaussian of roughly double the size of the original beam. With 
this smoothed image a mask was created including real emis- 
sion only. Thi s was checke d for all channels using the KARMA 
movie option (lGoochlI996h . This mask was used for CLEANing 
as described above for the self-calibration process of the contin- 
uum. This CLEANing procedure was repeated 3 times until no 
significant sidelobes were visible. Integrated intensity and ve- 
locity maps were created using the MOMENT task in MIRIAD. 
Finally, primary beam correction was applied to the integrated 
intensity maps. We created maps which were naturally weighted, 
thus providing the highest sensitivity, as well as maps with a ro- 
bust weighting of 0.5 to reach a higher spatial resolution. We 
reached an rms of 0.4 mJy beam"' in the naturally weighted 
map and and rms of 0.5 mJy beam"' in the robustly weighted 
map. Here we present only the robust map with the higher spa- 
tial resolution, because the naturally weighted maps, although 
more sensitive, do not show any additional significant flux con- 
tribution compared to the robust map. 

4. Results 

4. 1 . Neutral hydrogen 

4.1 .1 . Morphology of the tidal tails 

The neutral hydrogen distribution in NGC 4441 shows two 
prominent tidal tails to the north and to the south (see Fig.lTJ. 
The southern arm is blueshifted with respect to the systemic 
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Table 3. The H i properties of the total system, the tails and the 
central disc of NGC 4441. Heliocentric velocities are given. 



H I properties 


total system 






Vopt,hd (kms"') 




2674 


vel. range (km s"') 




2580 - 2900 


Fni (Jykm/s) 




4.83 


Mh,(10''Mo) 




1.46 


Mm/ Lb 




0.82 


central ring 






major axis (", kpc) 




32, 5.5 


minor axis (", kpc) 




19, 3.4 


Mhi(10**Mo) 




4.2 


northern tail 






vel. range (km s"') 




2700 ~ 2820 


extent (centre-tip) (' 


kpc) 


4.6, 48 


centre - clump (', kpc) 


1.5, 16 


Mh,(10«Mo) 




2.1 


southern tail 






vel. range (km s"') 




2580 - 2640 


extent (centre-tip) (' 


kpc) 


4.0, 42 


centre - clump (', kpc) 


1.7, 18 


Mh,(10«Mo) 




3.4 



velocity (see Fig|2]l, while the northern arm has a higher 
velocity. Both tails have a similar extent; the northern tail out to 
48 kpc, the southern tail out to 42 kpc from the nucleus and both 
form a relatively symmetric structure. The Hi tails are clumpy, 
including the tip of both tails. Two large clumps are found, both 
~ 17 kpc from the centre (marked as NC and SC in Fig.lTJ. 
From the kinematics there are no indications that these clumps 
are or will b e separate gravitationally bound entities like those 
simulated by lPuc & Bournaudl(l2007l) . 

It is noticeable that the symmetry in the Hi distribution is not 
seen in the optical morphology. Here, there is a tail to the north 
which coincides nicely in position and extension with the large 
northern H i clump. In the southern direction, however, we see 
a misalignment of the H i tail and the optical tidal features. Two 
shells are visible in the optical southwest of the main body. 
There is almost no overlap of the more western, outer shell with 
the H I tail. The tip of the inner shell corresponds to the southern 
Hi clump. Three optical knots are found at that position, of 
which at least two do not look point-like, i.e., are probably not 
foreground stars. They might be star forming knots, probably 
too close to the galaxy to survive as tidal dwarf galaxies, but can 
also be background galaxies, because our deep R-band image 
shows a background galaxy cluster However, no H i is detected 
at higher velocities within our band, which may of course miss 
distant galaxies. Unfortunately, the distance of the background 
cluster is not known. 

The gas in the tails shows only little dispersion in general 
(< 15 km s"'). The dispersion increases slightly in regions of 
gradients in the atomic gas distribution, e.g., around clumps. 
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Fig. 3. R-band image (taken with the CAHA 2.2 m telescope) of 
the central region. H i is overlayed as contours, contour levels 
are from 2.8 to 10.08 x lO^^'cm^^ in steps of 0.56 x lO^'^cm^^. 
The beamsize is shown in the upper left corner. The gas ring is 
centred on the optical nucleus, which exhibits a thin dust lane. 
The orientation of the dust lane is marked by the dashed line. In 
the lower left corner, the R-band image of the nucleus with the 
dust lane is shown. 

Table 4. Estimated parameters of the central Hi disc in 

NGC 4441. 



estimated parameters 




disc extent (") 


60 


max. velocity (km s"') 


163 


position angle (°) 


-35 


inclination (°) 


70 


Mdy„(10'"Mo) 


6.4 


Mm/Mj,.,, 


0.028 


Miyn/Ls (Mo/Lb,q) 


6.3 



concentration of H i is found on the west side of the ring, thus, 
west of the optical centre. Three H i peaks are present in the 
ring, with a distance of ~2.5kpc from the centre. This ring 
structure might actually be a small disc, where in the inner part 
most of the gas is present in the molecula r phase which is often 
found at column densities of ~ 10^' cm"^ (lYoungl2000l) . 

Fig. |4] shows the integrated Hi spectra of the total Hi amount 
and the central ring. The total Hi spectrum shows three peaks. 
The features around 2600 km s ' and 2780 km s ' correspond to 
the southern and northern tidal tail. A third peak at 2840 km s"' 
appears at the position of the northern clump (see also Fig.|2]i. 



4.1.2. The central ring 

Fig l3] shows the central region of NGC 4441 in the optical R 
band and the H i distribution (robust=0.5) overlayed as contours. 
We find a distinct H i ring with a diameter of ~ 1 1 kpc centred 
on the optical (and CO) nucleus. At the centre of the ring there 
is an H I depression, which is, however, unresolved. The highest 



4.1.3. Kinematics 

The main body of the H i distribution shows a relatively regu- 
lar rotation pattern which allows at least a crude estimate of the 
dynamical mass. Note that this region is much larger than the 
central ring seen in the H i distribution map and has a different 
position angle. 
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Fig. 4. Integrated H i spectrum of top: the total amount of H i in 
NGC4441 and bottom: the central ring. In the total Hi spec- 
trum, the two narrow features around 2600km s ' and 2780km 
s"' correspond to the two tidal tails. The third peak at 2840km 
s ' appears at the position of the northern clump. 



In Fig.|5] the region of the central velocity field is shown. 
Using the axis ratio it is possible to estimate geometrically the 
inclination. Even though drawing an ellipse is approximate, it is 
not arbitrary here. Based on the elliptical shape of the velocity 
field shown in Fig.|5] it is clear the the disc is seen relatively 
edge-on, thus the uncertainty of the inclination does not affect 
the mass calculation significantly. 

We measured an axis ratio of ~ 1/3, equivalent to an inclina- 
tion of ~ 70°. An error of +20° leads to an uncertainty of a factor 
of 2 in mass and thus gives at least an indication of the mass 
range. Note, that the minor axis is however not perpendicular to 
the major axis, which is a hint for non-radial motions like infall 
or outflow that could be triggered by the merger. One possible 
explanation would be an elongated gravitational potential, e.g., 
a bar or spiral arms, which forces the gas to move on ellipsoidal 
orbits. The resolution is, however, not sufficient to investigate 
this in more detail. Furthermore, the minor axis corresponds to 
the thin dust lane seen in the central region (see Fig. [3). 
In Fig.|6]we present the position-velocity-diagram (pv-diagram) 
through the major axis of the H i disc (position angle (pa) - 
325°). It shows a broad velocity range which is expected for 
a highly inclined disc. In the centre there is a gap because of 
the local minimum in the Hi distribution. A second velocity 
component at high velocities appears in the inner part, which 
demonstrates the complexity of the velocity distribution in this 
merger remnant. 
With an inclination of ~ 70° and a maximum Hi rotation 



64°50' 






L 


+49' 


- 


fl 


i 




- 




r 


+48' 


— 


— 


+47' 


- 


^ M^ 


- 




- 




- 


+46' 


II 


1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 



2700 ^ 



Jghg^m^QS 3Q5 2QS JqS 

Right Ascension (J2000) 



Fig. 5. Inner part of velocity field of NGC 4441, which shows a 
fairly regular solid body rotation pattern. The contour levels are 
from 2775 to 2825 km s"\ in steps of 25 km s"' 
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Fig. 6. Position-velocity diagram of the major axis of the inner 
H I disc in NGC 4441 . The position angle is 325°. Shown are the 
3 cr contours of 0.9 mJy beam"' . 



velocity of ~ 163 km s"' at a distance of r(estimated from the 
pv- diagram), we derive a dynamical mass of ~ 6 • 10'*' M© (see 
TableHJi. Even though this is a crude estimate, the mass-to-light 
ratio is within the range of what is found in normal galaxies. 
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Table 5. Continuum and FIR fluxes, SFR based on continuum 
and FIR emission, respectively, and q parameter. 



Right Ascension (J2000) 

Fig. 7. R-band image (CAHA) and contours of the 20 cm con- 
tinuum flux. The beam is shown in the lower left corner of the 
image. Contour levels are 0.5, 1, 2, 5, 7.5, lOmJybeam"'. 



4.2. 20 cm Continuum 

Fig. |7] shows an optical image with the 20 cm radio continuum 
overlaid as contours. The continuum source in NGC 4441 is 
barely resolved and concentrated on the nucleus of the galaxy. 
We fitted a Gaussian to the continuum source using the MIRIAD 
task IMFIT and derived a deconvolved source size of less than 
1 kpc. The derived centre coordinates corresponds to the optical 
coordinates. To estimate the star formation rate (SFR) from the 
20-cm radio continuum flux density we use 



SFR (Moyr-i) = 0.14 D^Faoc 



taken from lCondonl ( II992l) and iHaarsma et all ( |2000|) . where D 
is the distance in Mpc and F20 cm is the 20-cm radio continuum 
flux density in Jy. In Table |5]the results for NGC 4441 are given. 
Table |5]also gives the SFR estimated from the FIR luminosity; 

SFR (Moyr-i) = 0.17 Lfir 



following IkennicuttI (Il998l) . with LpiR in units of 10*^ Lp. The 
IRAS 60 Aim and lOO/im flux densities were taken from Moshitj 
(119901) . The SFR of 1-2 Mq yr ' calculated with both methods 
is low for a merger and normal for typical spirals. 
We also calculated the logarithmic ratio o f FIR to radio flux, the 
so-called q parameter (lHelouetalJll985h : 



? = log 



F(F1R)/(3.7510"Hz) 
F( 1.4 GHz) 



Note, that the continuum flux at 1.4 GHz, F( 1.4 GHz) 
in this formula, is in units of W m^ Hz ' - 
102''Jy. The FIR flux, F(FIR), is a combination of the 
60 /^m and 100/im flux densities, as used for determining 
Lfir : F(FIR) = 1.26 ■ lO^'^ . (2.58 ■ F(60yum) -1- F(100//m)). 
To calculate q, F(FIR) is normalised by the mean frequency of 
eOfim and 100 jum, which is 3.75 ■ 10'^ Hz. iHeTou et al.1 (Il985h 
derived a mean value for q of 2.3. Although this parameter was 
determined for disc galaxies only (but independent of their star 
formation activity), NGC 4441 (q = 2.4) conforms to this value. 



key characteristics 


20cm(Jy) 


0.013 


60;um (Jy) 


2.7 


lOOyum (Jy) 


3.9 


SFR2o,„,(M^yr') 


2.4 


SFi?FiR(Moyr-i) 


1.0 


1 


2.4 



5. Discussion 

5.1. The merger origin 

For an understanding of the whole merger process of galaxies 
resulting in a remnant like NGC 4441, it is necessary to know 
the progenitors. Of course it is difficult, if not impossible, to get 
a unique answer to that question from a completely merged rem- 
nant. However, there are some indications from our observations 
which at least narrow the parameter space. 
Since we found a large amount of gas in NGC 4441, at least one 
spiral galaxy must be involved in the merger process. In the fol- 
lowing, we discuss the possibility of a spiral H- spiral (S-hS) and 
spiral H-elliptical (Sh-E) merger 



Spiral+elliptical merger? 

There are several arguments that support an S+E merger origin: 

No strongly enhanced star formation: Typically, S+S merger 
remnants which show strong tidal features comparable to what 
is seen in NGC 4441, host a re gion of a (fading) starburst, e.g., 
Ikibbard & van Gorkoml (Il996l) . By looking at the FIR lumi- 
nosity, the 20 cm con tinuum emission and the optical spectra 
(iManthev et al.l 120051) . however we do not find any indication 
for a strong starburst, neither ongoing nor in the recent history, 
in NGC 4441. Therefore, less gas is probably involved in this 
merger event, which is naturally the case when one progenitor 
is an elliptical. 

It has to be mentioned, however, that the strength of an induced 
starburst not only depends on the available amount of gas 
but also on the internal structure of the progenitor galaxies. 
In particular, the presence of a bulge in a progenitor spiral 
galaxy will stabilise th e gas disc aga inst the formation of a 
bar during the merger (lMihosl[l999al) . That bar would be a 
powerful mechanism to support in fall of gas into the centra l 
region leading to a strong starburst (Miho s & Hernquisall994l) . 
In S-hS merger simulations gas dissipates energy and thus 
loses a ngular momentum when it becomes shocked along 
a bar ("Barnes & HernquistI ll99lL 119961: iMihos & Hernquisli 
[l996; Barnes 1993)- This leads to rapid gas inflows and a 
concentration of gas in the centre of the remnant, resulting in a 
starburst. 

Nevertheless, the simulations for S-hS with bulges, i.e., no 
merger-induced bars, never show such a low star formation rate 
as in the case of NGC 4441 . 

Of course, it is necessary to investigate in numerical simulations 
how bar formation and a present bulge influences the gas 
in a S+E merger scenario, in particular when the galaxies 
are completely merged. Due to the spheroidal shape of the 
elliptical which is similar to a bulge, bar formation is likely not 
important and therefore no strong gas inflows occur Since the 
major and minor axes of NGC 4441 are not perpendicular, this 
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might indicate the presence of a small oval distortion which is, 
however, not strong enough to trigger enhanced star formation. 

Theoretical models: Various numerical simulations showed 
that optical shells in combi nation with ti dal tails are typical 
for S + E merger remnants ( Ouinn 19841: iD upraz & Combes 
1986 f: 'Hernquist & Ouinn' 1988; 'Kojima & Noguchi| I1997t 
Combes & Charmandaris 1999). These shells are created by 



phase wrapping of th e disc in the potential well of the elliptical 
galaxy (lOuinnll 19841) . They occur in the advanced merger stage 
~ 4-5 ■ 10^ yr after the impact and last for more than 1 Gyr 
dQuinh 1984; Kojima & Noguchi 1997D. 

The simulated gas distribution stron gly depends on the model 
used f or th e ISM. It has been noted bylCombes & CharmandarisI 
( 119991) and lCharmandaris & CombesI (1200 0^ that the often used 
smoothed particle hydrodynamics (SPH) approach might not 
always be appropriate. Instead, these authors account for the 
clumpy structure of the ISM by simulating the behaviour of 
individual gas clouds rather than a continuous fluid. They found 
a largely spread gas distribution, often forming gaseous tails as 
we see in NGC 4441. Therefore, these models do not predict a 
strong starburst phase, because the gas becomes too dispersed. 
This agrees with our obser vations, but stands in con t rast to 
S+S m erger scenarios (e.g., iBarnes & Hernquislll991l; iMihosI 
[HUB). 

Stellar component similar as the Medusa: The similarity of 
flie optical morphologies of NGC 4441 and NGC 4194, the 
prototypical S+E merger remnant, has been mentioned in the 
introduction, reflecting a possible analogue merger origin. 
Besides the obvious features like the stellar shells and tidal tails, 
both galaxies host a dust lane. 

In NGC 4441, the minor axis of the fairly regularly rotating 
Hi disc coiTesponds to a faint central dust lane (Fig.[3]l. It is 
interesting to note that also the Medusa exhibits a minor axis 
dust lane. Although this may further strengthen the notion of 
a a similar origin, the dust lane of the Medusa is more promi- 
nent and longer (it curves into the optical tail). Furthermore, 
t hat dust lane is associat ed with a double spectral feature 
JAalto & Htittemeisteij 12000 ). Such extra features are not seen 
in NGC 4441. This may be a result of its more evolved state, 
because a large fraction of the ISM has been consumed in star 
formation. 

Hi vs. optical morphology: In NGC 4441 we see striking differ- 
ences between the stellar and gaseous distribution. In particular 
the two optical shells stand in contrast to the two symmetric H i 
tails. 

In S+S merger remnants, h owever, the optical and H i mor phol- 
ogy are generally similar dHibbard & van GorkomI Il996h . be- 
cause both the stellar and gaseous tidal fe atures ar e formed from 
the two discs which are merging. Hibba rdet al.l yOOO) found 
some examples in which the optical and the gaseous tails are 
displaced or where some regions are lacking gas, while the over- 
all morphologies are still similar. The authors described several 
scenarios for these discrepancies between the stellar and gaseous 
tidal morphology and claim that mainly superwinds are respon- 
sible for the differences found. These superwinds are driven by 
the central starburst occuring in S+S merger remnants and sweep 
away the gas in the direction of the wind. 
Their findings differ, however, in two ways from NGC 4441. 
First, NGC 4441 does not host an ongoing strong starburst, nei- 
ther did it in the recent history. Therefore, the existence of a 
superwind is unlikely. Second, we see not only an offset of the 



gaseous and stellar tidal features, but also a completely differ- 
ent morphology. This cannot be explained by a possible super- 
wind alone. Other possible explanations given by Hibbard et al] 
(l2000h . like dust-obscuration and photoionisation, can only ex- 
plain local differences and not large-scale discrepancies, as 
pointed out by the authors. Photoionisation might explain why 
there is only little gas found in the outer optical shell. In that 
case we expect to find a large amount of ionised gas which can 
be detected in e.g. Ha emission. This will be done in a further 
investigation. 



Spiral+Spiral merger? 

There is one striking difference between NGC 4441 and S+E 
merger simulations as well as the prototype NGC 4194: 

Two symmetric Hi tails: The discussion of the merger origin 
of N GC4441 is complicated by the finding of two sym- 
metric Hi tails. This is typical for a merger between two 
disc galaxies, when each disc forms on e tail during the en- 



counter (e.g., Tbomre & Toomrel 119721: iHernquist & Barnes! 



1994; Mihos & Hernguisll 1 19961) . Neither the sini u lations 



by iWeil & Hernguisll (1 19931) nor |Kojima & Noguchi| (1 19971) 
show the appearance of two tails in an S+E merger scenario. 
Furthermore, our Hi data of the prototypical S+E merger 
NGC 4 194 show only one gaseous tail in that galaxy, in accor- 
dance to the S+E merger simulations (Manthey et al. 2008b) 
Observations of the interacting S+E galaxy system Arp 140 
(an early stage of an S+E m erger where both g alaxies are 
still clearly separated) done by ICuUen et alj ( l2007h showed as 
well one Hi tail only. Furthermore, they found an Hi bridge 
between the interacting galaxies formed by disc material 
which streams towards the elliptical galaxy dToomre & Toomrg 
1972) . This bridge will likely be destroyed once the galaxies 
finally merge, as in the case of bridges found in S+S mergers 



In rece nt deep observations of early-type galaxies lMorganti et al.l 
("2006") found a significant amount of H i in a large fraction of 
those galaxies. The data reveal a variety of structures, includ- 
ing extended discs. However, the peak column density is with 
N(H i) < 10^" cm"^ lower than what we found in NGC 4441 
(N(H i) - 1.1 ■ 10^' cm"^), so no sta r formation oc c urs in 
the regular early-type galaxies of the iMorganti et al.l d2006h 
sample. However, considering the likely S+E merger origin of 
NGC 4441, it is possible that the elliptical progenitor possessed 
such a gas disc, which would explain why NGC 4441 has two 
H I tidal tails and otherwise resembles an S+E merger remnant. 
The low column density Hi disc around ellipticals might also 
explain why the star formation in NGC 4441 is not as enhanced 
by the merger as in 'normal' disc-disc mergers. 

5.2. The formation of the central Hi disc 

Clearly, NGC 4441 is in an advanced merger stage, because 
even the progenitor nuclei are completely merged in addition to 
the pronounced stellar tidal features. In this stage a stable gas 
disc can already have formed, as shown by iBarnesI d2002h . In 
unequal-mass mergers the disc formation is more impressive 
dBarnesll20^ . which again supports that NGC 4441 is at least 
not an advanced remnant of a merger between two large disc 
galaxies. 
The fairly regular rotation pattern in the central HI disc indicates 
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that the merge r event has h appened several hundred milhon 
years ago (e.g..' Mihoslll999a) . The gas in this advanced merger 
has already started to settle down, but without reaching an 
equilibrium state yet. The lack of Hi in the centre (Figs.[Tl |6]l 
might be caused by an efficient transformation of atomic to 
molecular gas . This is supported b y our CO detection within the 
central 22"(Mant hevet 



jpported p ' 
ai] l2008al) ." 



The velocity gradient of the central disc connects to that of the 
two tidal tails, even though perturbations are clearly present. 
That means that the material in the tails and the central disc 
share the same sense of rotation. This indicates that a large 
amount of angular momentum still remains in the gas of the 
tidal tails. A merger geometry in which the angular momentum 
does not cancel out much seems therefore likely, e.g., a prograde 
merger with parallel rotation spins of the galaxies. Independent 
of the presence of an Hi disc in the elliptical progenitor, a 
prograde merger is likely in which the rotation of the spiral and 
the movement of the spiral towards the elliptical are aligned. 
The amount of conserved angular momentum in the gas also 
depends on other factors, e.g., the presence of a bar formed 
during the interaction in which the gas loses angular momentum 
by shocks (see above). 

If no strong shocks occur, the gas retains much of its angular 
momentum, which prevents the gas from falling toward the 
central region on short timescales. Instead, a gaseous disc will 
be formed. The remaining gas in the tidal tails will finally 
return to the main body and will lead to the growth of the 
present Hi disc. This behaviour is modelled by e.g., iBarnesI 
( 2OO2I) and is observed in other merger remn ants like NGC 7252 
(Hibbard et al.'l994' 'H ibbard & MihosI 19950 and NGC 4038/39 
( Hibbardetal. 2001 ). 



5.3. The star forming region 

NGC 4441 has a moderate ongoing star formation rate of 
1 - 2Mo yr '. However, there are indications based on optical 
spectroscopy that ther e has been a period of enhanced star 
formation ~ 1 Gyr ago ( lBergvallll98ll:lManthev et al.ll2005h . 
The extent of the star-forming region can be investigated 
using the continuum data. The deconvolved size of the 20 cm 
continuum source was determined as <lkpc (Table|5]l. The 
emission comes likely from ongoing star formation only, since 
the q value agrees with the FIR-to-20cm continuum ratio 
typical for star formation, and based on optical spe ctroscopy 
there a re no hints for another origin, e.g. an AGN (iBergvalll 
Il981h . iManthev et al.l (l2008al) found molecular gas extended 
out to ~4kpc in NGC 4441. Thus, the region of ongoing star 
formation is embedded in the reservoir of raw material for 
star formation, namely the molecular gas extent, but does not 
cover its full extent. Compared to the Medusa, where the region 
of ongoing star formation is extended on a scale of ~ 2kpc 
(e.g.. lArmus et al. 1990) and much more intense (SFRopt ~ 
30 M0 yr~^), the star forming activity seems to be shrunk in 
intensity and maybe also in spatial extent. The question arises, 
why the starburst in NGC 4441 has alrea dy faded, while the 
Medusa is still intensely forming stars. Aalto & Hiittemeisteij 
(I2000.) estimate that the burst in the Medusa will continue only 
for ~ 40-10^ years, under the assumption that the current 
star formation rate will not change dramatically. If the star 
formation rate depends mainly on the duration of the burst 
and decreases with ongoing burst age, the phase of intense star 
formation in a merger event, which NGC 4441 and the Medusa 



both underwent, might be very short, and the Medusa will fade 
soon as well. 

Furthermore, it has become clear by simulations of S+S merg- 
ers that the starburst intensity in mergers of unequal masses 
strongly depends on the structure of the progenitors (e.g., the 
presence of a bulge) and the geometry of the merger even t 
(orbits, spin orientation of the galaxies etc.) dCox et al.ll2007h . 
Even though no detailed simulations exist to investigate this 
phenomenon in the case of Sh-E mergers, it is likely that these 
factors play a similar important role in S+E mergers. Therefore, 
(unequal-mass) S+E mergers are expected to show a much 
larger scatter in starburst strength even if they are in the same 
merger stage. 



6. Summary 

1. We observed the moderate luminosity merger NGC 4441 
with the Westerbork radio interferometer and found a com- 
plex, extended H i structure affected by tidal forces. The total 
Hi mass adds up to 1.46 ■ 10^ Mq. 

2. We detected two Hi tails. Both tails have a similar extent 
out to 4'(42 kpc, southern tail) and 4.6'(48 kpc, northern tail) 
away from the centre and they emerge from opposite sides 
of the main galactic body. 

3. The northern H i arm follows the optical tail closely, however 
the H I is more extended and shifted to the east. In contrast, 
the southern H i arm is not aligned with the shell structure 
found in the optical, but touches only the inner optical shell, 
whereas in the outer shell no neutral gas is found. 

4. The velocity field shows a fairly regular rotation pattern in 
the inner region, associated with the optical main body. The 
overall velocity gradient is continued in the tidal tails. The 
southern tail has a generally lower velocity as the system, the 
velocity of the northern tail is larger than the central velocity. 

5. From the 20 cm continuum flux we calculated the star for- 
mation rate as 2.4 M© yr"'. From the FIR fluxes we get a 
SFR of 1.0 M0 yr"'. The resulting ratio q is 2.4, which is in 
good agreement with the typical q value in the literature. We 
determined a continuum source size of less than 1 kpc. 

6. Because of the large differences in the stellar and gaseous 
distribution, it is unlikely that NGC 4441 is a merger be- 
tween two large disc galaxies. If diff'erences are found in 
those mergers, they are on much s maller scales and can be 
best explained by superwinds (see iHibbard et al.ll2000h . In 
NGC 4441, no such strong galactic winds are expected, be- 
cause the galaxy is in a phase of moderate star formation. A 
merger between two ellipticals can be also excluded, because 
we found a large amount of atomic and molecular gas and 
signs of a past strong star formation activity. This properties 
are not expected in E+E mergers. Thus, it is most likely, that 
NGC 4441 is the remnant of an S+E merger. 
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Fig. 1. NGC 4441 - top left: H i distribution overlayed on optical DSS image. Contour levels are 2.8, 8.4, 14, 28, 56, 84 x 10''' cm"^ 
Marked are the optical features: the optical tail (OT), the inner shell (IS), the outer shell (OS) and the three optical knots, top right 
H I distribution, contours same as above. Marked are the gaseous features as described in the text: the northern tail (NT), the northern 
clump (NC), the central ring (CR), the southern clump (SC), and the southern tail (ST), bottom left: velocity field, contour levels 
are 2570 to 2840km s"', in steps of 30km s"', bottom right: Hi velocity dispersion (2nd moment). The cross marks the optical 
centre position. 
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Fig. 2. Channel maps (robust weighted) of NGC 4441. The optical centre position is marked with a cross. The contour levels are 
-0.9,1.1, 1.3, 1.5, 1.7, 1.8, 2mJy beam"'. The heliocentric velocity is given at the top left, the beam at the bottom left corner. 



